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ABSTRACT 

Numerical simulations of dwarf spheroidal galaxies undergoing sev¬ 
eral close encounters with the Milky Way are described. By compar¬ 
ing our models to observed properties of the recently discovered dwarf 
galaxy in Sagittarius (Sgr), we discuss implications of our results for 
the formation and evolution of the Milky Way system. We find that 
existing observations are not sufficient to allow us to place precise lim¬ 
its on either the orbit or the initial state of the dwarf. Debris from the 
ongoing tidal stripping of the Sagittarius galaxy are expected to form 
moving groups in the halo of the Galaxy and the discovery of such stars 
would strongly constrain the history and dynamical state of the dwarf. 
Furthermore, if Sgr is presently being disrupted, we predict that its 
remains will be detectable as a moving group in the halo for more than 
1 Gyr. Thus, if similar accretion events have occurred in the recent 
history of the Galaxy, their aftereffects may still be observable. 
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1. Introduction 

In April 1994, Ibatha, Gilmore & Irwin 
(1994; hereafter IGI) announced the discovery 
of a moving group in the direction of Sagit¬ 
tarius, centered on the Galactic coordinates 
1= 7.5°, b= —15°, and which is 24 ± 2 kpc 
from the Sun (~ 16 kpc from the center of the 
Galaxy). Surface density contours of this ob¬ 
ject are elongated perpendicular to the disk 
of the Milky Way and are characterized by 
an axis ratio of ~ 3 : 1. The lowest and 
highest contour levels, which cover the region 
/ = 5° to 10°, b = —20° to —10°, span a fac¬ 
tor of ~ 10 in surface density. The stars in 
the group have a heliocentric radial velocity 
of 153 ± 2 km/s, corresponding to 176 km/s 
in Galactocentric coordinates, a velocity gra¬ 
dient of < 1 km/s per degree along the major 
axis, and a velocity dispersion of 10 km/s. 

Mateo et al. (1994) subsequently obtained 
deep photometry of a held in this area and es¬ 
timate that members of the group have typ¬ 
ical ages of ~ 10 Gyrs and metallicities of 
[Fe/H]ss —1.0 ± 0.3. These results, together 
with the lack of HI in this region (IGI), sup¬ 
port IGI’s conjecture that the moving group 
is a dwarf spheroidal galaxy (dSph), which 
would make it the ninth thus far discovered 
around the Milky Way. IGI compare Sgr di¬ 
rectly to the Fornax dSph since the number 
of horizontal branch (HB) stars they detect 
(~ 17000) suggest that the two dwarfs are 
similar in size, leading to a magnitude esti¬ 
mate of M v ~ —13 (or L tot ~ 10 7 T@) for Sgr. 

Sgr’s current location offers a rare oppor¬ 
tunity to study the effects of a strong tidal 
interaction in detail. The tidal radii of clus¬ 
ters of stars of masses 10 8 M@ and 10 7 M@ at 
16 kpc from the center of the Galaxy are 0.77 
kpc and 0.36 kpc, respectively, assuming an 
enclosed mass of 3 X 10 n M@, whereas the 


group extends over ~ 2 kpc, suggesting that 
the dwarf is currently being disrupted. The 
extreme axis ratios inferred for Sgr are a re¬ 
flection not only of the compressive shock of 
the encounter on the dwarf but also of the 
dispersal of stripped material along its orbit. 
Thus, the direction of Sgr’s elongation traces 
its projected orbit. 

We model the dynamics of the Milky Way 
and Sgr dwarf by considering encounters that 
are consistent with the observed properties of 
the group. We examine whether or not ex¬ 
isting observations are sufficient to determine 
the dynamical state of the group; namely, if 
it is bound or unbound and whether or not it 
has suffered previous close encounters. Since 
only one component of the dwarf’s velocity is 
known, we do not know its precise orbit and 
cannot hope to reproduce all observations ex¬ 
actly, but rather must be content to inves¬ 
tigate the general characteristics of such en¬ 
counters. 

The mere existence of Sgr has implica¬ 
tions for models of galaxy formation since the 
timescale for the encounter is ~ 16 kpc/176 
km/s < 10 8 years. Clearly, accretions of 
dwarf galaxies by the Milky Way must be 
common or we must instead conclude that 
we are observing the Galaxy during a special 
phase of its evolution. The former possibil¬ 
ity supports Searle and Zinn’s (1978) proposal 
that the Galaxy was slowly built up from pri¬ 
mordial fragments and implies that the known 
dSph’s represent only a small fraction of a 
much larger original population of satellites 
(for reviews see Larson [1990] and Majewski 
[1993a]). To investigate this possibility, we 
follow the evolution of one simulation beyond 
complete tidal disruption and discuss observ¬ 
able consequences of such an event. 
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2. Method 
2.1. Models 


A three-component model is used for the 
Galaxy, in which the disk is represented by a 
Miyamoto-Nagai potential (1975), the spheroid 
by a Hernquist potential (1990a) and the halo 
by a logarithmic potential: 
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Here, M disk = 1.0 x 10 11 , M spher = 3.4 x 
lO 10 , Vhaio = 212, a = 6.5, b = 0.26, c = 0.7, 
and d = 12.0, where masses are in M®, ve¬ 
locities are in km/s and lengths are in kpc. 
This choice of parameters provides a nearly 
flat rotation curve between 1 and 30 kpc and 
a disk scale height of 0.2 kpc. The radial de¬ 
pendence of the z epicyclic frequency (n z ) in 
the disk between radii at 3 and 20 kpc is sim¬ 
ilar to that of an exponential disk with a 4 
kpc scale length. 

Initially, the satellite is represented by a 
Plummer model 
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summarizes each of the simulations and pro¬ 
vides analytic estimates of tidal radii (King 
1962) and the impulsive energy input for each 
model on its first encounter (see Binney & 
Tremaine 1987). Note that Model D em¬ 
ployed 50 times as many particles as the oth¬ 
ers to permit a more detailed analysis of phase 
space. 

Three different orbits were used in the sim¬ 
ulations. Table 2 summarizes the properties 
of these orbits and Figure 1 traces their paths. 
Each orbit is constrained to have a Galacto- 
centric radial velocity of 176 km/s at a dis¬ 
tance of 16 kpc from the center of the Galaxy 
and is chosen to lie in a plane perpendicular 
to the disk, in view of Sgr’s direction of elon¬ 
gation. Orbits 1 and 3 illustrate the effect of 
varying eccentricity, assuming that Sgr has al¬ 
ways been bound to the Milky Way, and are 
run in our full Galactic potential. An orbit 
much less eccentric than Orbit 1 would lead to 
the complete destruction of our dwarfs in un¬ 
acceptably short (< lGyr) timescales, while 
an orbit much more eccentric than Orbit 3 
would have an apocenter beyond the most dis¬ 
tant observed dSphs (Leo I and Leo II). Orbit 
2 is virtually identical to Orbit 1, but includes 
the effect of only the spheroid and halo poten¬ 
tials so that the relative influence of the disk 
can be determined. 

2.2. Evolution Code 


where r 0 = 0.6 kpc and the mass is either 
M = ll) 7 .\/. or 4/ ll)\\/.. 

For r 0 = 0.6 kpc, a mass of 10 7 M@ yields a 
central density comparable to the luminous 
component of the Fornax dSph and a cen¬ 
tral velocity dispersion of 3.4 km/s, while the 
higher mass includes the full dynamical mass 
of Fornax (Mateo et al. 1993) and gives a cen¬ 
tral velocity dispersion of 10.7 km/s. Table 1 


All computations were performed on the 
Connection Machine 5 (CM-5) at the Na¬ 
tional Center for Supercomputing Applica¬ 
tions (NCSA). The Milky Way is represented 
by a rigid potential, while the dwarf is de¬ 
scribed by a collection of self-gravitating par¬ 
ticles whose mutual interactions are calcu¬ 
lated using a self-consistent held code (Hern¬ 
quist & Ostriker 1992). Our code has been 
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parallelized to run on the CM-5 (Hernquist, 
Sigurdsson & Bryan 1994) to maximize par¬ 
ticle number and consequently minimize nu¬ 
merical relaxation (Hernquist & Barnes 1990, 
Weinberg 1993, Johnston & Hernquist 1994). 
Since the mass of the Sgr dwarf is much 
smaller than the mass of the Milky Way, dy¬ 
namical friction and energy exchange between 
the parent and satellite galaxies are assumed 
negligible and have been ignored. 

2.3. Analysis 

To interpret our results, we convert the 
contours of the observed isopleth map of Sgr 
(see IGI) into mass surface densities. The 
lowest contour seen is at 0.125 image/(arcmin) 2 , 
which corresponds to 0.0026 image/pc 2 at a 
distance of 24 kpc, where an “image” is an HB 
star. This number density, n, can be related 
to a mass surface density using 



where A is the total luminosity of the system 
per red giant branch (RGB) star with magni¬ 
tudes less then the magnitude of the HB, B 
is the ratio of the number of HB stars to the 
number of RGB stars with magnitudes less 
than the magnitude of the HB, and C is the 
mass-to-light ratio for the system. The mag¬ 
nitude of the HB can be approximated by 

M hb = 0.66 + 0A9[Fe/H] 

(Bolte 1994), and A can be calculated from an 
appropriate luminosity function. Using Berg- 
busch and VandenBerg’s (1992) isochrones for 
an 8.0Gyr, [Fe/H] = —1.03 cluster, and an 
initial mass function exponent x = 1.5 gives 
A = 528. B can be estimated from observed 
luminosity functions and ranges from B = 1 
to B = 2 for various ages and metallicities 


(e.g. Mighell 1990; Buonanno, Corsi & Fusi 
Pesci 1990). C is 2.25 for the chosen lu¬ 
minosity function. Hence we estimate that 
E > lM@/pc 2 for the lowest contour in IGI’s 
isopleth map. 

3. Results 

In what follows we first consider the ef¬ 
fects of varying model parameters and then 
compare results to observations of Sgr. Fig¬ 
ure 2, described in §3.1, illustrates the evolu¬ 
tion of each model; Figures 3-7, described in 
§3.2 — §3.4, summarize the state of the models 
at times consistent with Sgr’s current position 
and velocity relative to the Sun; and Figures 
8 - 10, described in §3.5 — §3.7, explore the 
existence and persistence of moving groups in 
the halo. 

Observations of Sgr constrain the average 
Galactocentric radius (16 kpc) Galactocentric 
radial velocity (176 km/s), and position rel¬ 
ative to the Galactic disk (15°) of the most 
tightly bound particles in the simulation. To 
exactly match all these parameters on each of 
many passages along an orbit with a given 
pericenter and apocenter would require a cor¬ 
responding number of simulations, with each 
one chosen to match the constraints on one 
particular passage, because the orbit is a 
rosette and hence the angular position of peri¬ 
center precesses. However, the halo and the 
disk contribute nearly equally to the disrup¬ 
tion (see §3.1) and the dwarf receives roughly 
the same energy impulse on each encounter 
(to within a factor of two) irrespective of the 
exact orientation of the disk relative to the 
orbit or location of disk passage. Indeed, Fig¬ 
ure 6, which illustrates the fifth passage of 
Model D, shows no anomalous features when 
compared to the other figures, even though 
its most recent disk passage occurred at more 
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than 30 kpc and the impulse from the disk was 
significantly smaller than in the other pas¬ 
sages. In addition, the current position, ra¬ 
dial velocity and elongation of Sgr itself sug¬ 
gest that it is not far from its own pericenter 
and disk passage and is likely to have suffered 
from a similar impulse to that experienced by 
the models in our simulations. Hence, as a 
practical compromise, we drop the last ob¬ 
servational constraint (orientation relative to 
the plane of the Galaxy) and, in Figures 3-7, 
consider the state of the dwarf only when it is 
at a radius of ~ 16 kpc with a radial velocity 
of ~ 176 km/s. 

3.1. Mass Loss and Disruption 

Figure 2 shows the fractional mass still 
bound by self-gravity as a function of time 
during each of the simulations. The curves 
clearly show that mass loss occurs primarily 
and almost instantaneously at pericenter ow¬ 
ing to the impulsive energy input from the 
encounter. The orientation of the disk with 
respect to the orbit and the exact distance 
of the dwarf from the Galactic center as it 
passes through the disk, which both vary on 
each passage, do not affect the impulse sig¬ 
nificantly since the mass lost at pericenter is 
always of the same order. 

By comparing Models B and D, which dif¬ 
fer only in their orbits, we see that the mass 
lost at pericenter depends weakly on orbit pa¬ 
rameters. Model B survives fewer encounters 
than Model D, however, because it loses mass 
at a greater rate throughout its (less eccen¬ 
tric) orbit. 

Models B and C, which differ only in that 
the Galactic disk is excluded in Model C, 
show substantial differences in the rate of 
mass loss. This result emphasizes that rather 
than being merely a disk-shocking phenomenon, 


tidal disruption involves both disk and halo 
shocking, the latter being much like bulge- 
shocking of globular clusters (e.g. Aguilar, 
Hut & Ostriker 1988). In fact, the internal 
potential energy of bound particles in Model 
C after its fourth encounter is 1/3 its initial 
value, representing an energy gain of 2/3 its 
initial binding energy. Since Model B, run in 
the full potential, was completely disrupted 
over the same number of encounters, we infer 
that the last 1/3 of its binding energy must 
have been gained from the disk. We conclude 
that at least half of the energy gained from the 
encounters in the full potential comes from 
halo shocking alone. In addition, from the 
distorted shape of its surface density contours 
it is impossible to infer if Sgr has yet passed 
through the disk since the disk is not solely 
responsible for the satellite’s destruction. 

3.2. Surface Densities 

In Figures 3-7 we present results from the 
viewpoint of an observer at rest with respect 
to the center of the Galaxy and 24 kpc from 
the dwarf on the opposite side of the bulge. 

Panel (a) of each figure shows mass sur¬ 
face density contours, with the bold contour 
at the level lM@/pc 2 . The isopleths of the ob¬ 
served star counts are much less smooth than 
the contours in our figures, probably due to 
uncertainties in the counts themselves and be¬ 
cause of dust obscuration towards the bulge 
of the Galaxy. However, each simulated map 
compares well in extent, ellipticity, density 
contrast, and absolute density to the obser¬ 
vations. 

Panel (f) shows the surface density profile 
calculated on concentric ellipses, with an el¬ 
lipticity chosen to fit the contour at 5° along 
the major axis (solid line). The dotted line is 
a King model of the given concentration, plot- 
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ted for comparison. The central surface den¬ 
sity and concentration parameter of the King 
model fitted to the profiles both decrease with 
increasing number of pericentric passages un¬ 
til a critical point is reached, after which full 
disruption occurs. On the last passage prior 
to disruption, all our models (irrespective of 
orbit and initial conditions) are well fitted by 
nearly identical King models. This interesting 
coincidence merits further investigation, but 
is not something we discuss further here. A 
fuller exploration of parameter space (varying 
both the mass and the scale length of the ini¬ 
tial model) would be necessary to address this 
issue since the criterion for tidal disruption 
involves the density contrast of the satellite 
and parent galaxies. (Previous numerical in¬ 
vestigations of tidal interactions and density 
criteria for disruption include Miller [1986], 
McGlynn [1990], and Oh, Lin and Aarseth 
[1994].) 

3.3. Kinematics 

For the kinematic analyses all velocities 
are projected along the line of sight noted 
in the previous section. Particles are or¬ 
dered by their position along the major axis 
(positive angles are in the direction of mo¬ 
tion of the dwarf) and then binned into eight 
(1.25°) 2 boxes before calculating the velocity 
moments. No bin contains fewer than 100 
particles. 

Panel (b) shows the average velocity (open 
symbols) and velocity dispersion (filled sym¬ 
bols) of particles in bins along the major axis 
of the contour map. ft is interesting to note 
that tidal stripping does not enhance the line- 
of-sight velocity dispersion, even out to pro¬ 
jected distances of 10°. 

The models do have a substantial positive 
velocity gradient (> 2 km/s per degree) along 


their major axes, in contrast to the observa¬ 
tions of < 1 km/s per degree. Two effects 
contribute to line of sight velocity gradients 
observed along the projected major axis of a 
satellite in the tidal held of the Milky Way. 
The first is due to the combination of the 
tangential velocity, v t} of the satellite with 
its angular extent on the sky. An angular 
separation of A 6 between two lines of sight 
along the major axis leads to a velocity dif¬ 
ference of rysinA# (assuming that the major 
axis is coincident with the tangential veloc¬ 
ity vector). The second is due to stripping 
of material from the satellite. Stars form¬ 
ing the leading streamer are pulled towards 
the Galaxy to lower energy orbits and those 
forming the trailing streamer are pulled away 
to higher energy orbits (Oh, Lin, & Aarseth 
1994). Clearly the sign and magnitude of the 
gradient will depend on the position of the 
observer relative to the Galactocentric radius 
vector of the satellite and the strength of the 
tidal interaction. In our case (for a line of 
sight coincident with the radius vector), the 
first effect leads to a positive velocity gradi¬ 
ent in the direction of motion of the satellite, 
while the second leads to a negative gradi¬ 
ent. In Figures 3-7 the dwarf has a large 
angular extent on the sky and the first effect 
dominates, with Model D (Figures 6 and 7) 
showing the most significant gradient because 
it is on a higher energy orbit and has a larger 
v t . That such a gradient is not observed in 
Sgr could be because the two effects are of 
the same magnitude. 

Panel (c) shows normalized values of the 
Gauss-Hermite moments calculated for par¬ 
ticles in each bin. These moments quantify 
the deviation of the velocity distribution from 
a Gaussian (Gerhard 1993; van der Marel & 
Franx 1993). We have chosen van der Marel 
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& Franx’s formalism using the derivation out¬ 
lined in Heyl, Hernquist & Spergel (1994). 
Since the higher order moments are only a 
few percent of the zeroth order term, the dis¬ 
tribution is nearly Gaussian along the major 
axis, as confirmed in panels (d) and (e) which 
show the velocity distribution of particles in 
the bins at —9° and —1°, respectively. 

3.4. Comparison with Observations 

There are no striking differences between 
the results shown in Figures 3-7, and each 
fits most, but not all, of the observations well. 
The points of ambiguity are the smoothness 
of the density contours and the large veloc¬ 
ity gradients seen along the projected major 
axes in our simulations. The absence of such 
a velocity gradient in the observations is very 
surprising in the light of these simulations and 
may provide an important clue to Sgr’s his¬ 
tory. If the lack of gradient is confirmed in 
further observations then it will place strong 
constraints on the orbit of Sgr by requiring 
that the two effects contributing to any ve¬ 
locity difference (as described in §3.3) exactly 
cancel. 

In summary, our results show that Sgr is 
neither necessarily on its first nor on its last 
pericentric passage. There are no outstanding 
features in the velocity distribution to distin¬ 
guish between a passage in which a bound 
core of particles still remains, as in Figures 
4 and 6, from one in which all particles are 
unbound, as in Figures 5 and 7. Severely dis¬ 
torted contours are apparent in all these fig¬ 
ures. Thus it is difficult to place constraints 
on the initial mass of the dwarf or on its orbit 
with current observational data. 


3.5. Phase Space Structure 

For the remainder of §3 we discuss the 
position-velocity phase space evolution of Model 
D in detail. This model employed 50 times as 
many particles as the others to improve res¬ 
olution and accuracy. In addition, this was 
the only model with a disruption timescale of 
order the Hubble time. 

Figure 8 shows the radial velocity as a func¬ 
tion of radius for a sample of 2500 particles 
immediately following the final passage of the 
satellite (at the time illustrated in Figure 7). 
The solid line traces the dwarf’s orbit. At 
this time the model has undergone 6 pericen¬ 
tric passages over the course of 13 Gyrs and 
it has been fully disrupted. Each encounter 
has produced its own distinct tidal streamer. 
Nevertheless, the particles occupy a relatively 
small region of phase space. We hope to ex¬ 
ploit this high density to locate tidal debris 
from cannibalized dwarfs in the halo of our 
Galaxy. 

3.6. Streamers as Moving Groups 

The Sagittarius dwarf was first recognized 
as a moving group. Is it possible to detect 
the streamers from the encounter away from 
the Galactic center? In Figure 9, we exam¬ 
ine the line-of-sight velocity distributions of 
stars in two different square-degree regions. 
The viewpoint is again that of an observer at 
rest with respect to the center of the Galaxy 
and 24 kpc from the dwarf on the opposite 
side of the bulge. Panels (a) and (b) are for 
stars 40° away from the dwarf along its or¬ 
bit, while panels (c) and (d) show similar in¬ 
formation at 75°. Apparent magnitudes are 
estimated using the adopted luminosity func¬ 
tion (solid line). Also plotted in each of the 
panels is an estimate of the observed numbers 
of stars as a function of line-of-sight veloc- 
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ity in these directions (dotted line). These 
distributions were generated using the star 
counts of thin disk, thick disk, and halo stars 
in the direction of the North Galactic Pole 
(Reid & Majewski 1993). The velocity dis¬ 
tribution of each component is modeled as a 
Maxwellian with z-velocity dispersions of 30, 
45, and 150 km/s for the three components, 
respectively. A comparison of the areas un¬ 
der the two curves in each panel shows that 
the surface density of the stars in the stream¬ 
ers is only ~ 10% of the background. How¬ 
ever, the tidal streamers stand out as mov¬ 
ing groups in all magnitude bins, though the 
exact distribution depends on the details of 
the encounter. Observing such a group di¬ 
rectly associated with the Sagittarius dwarf 
would strongly constrain its orbit, and hence 
its history. The tidal streamers are expected 
to trace the great circle defined by Sgr’s elon¬ 
gation. 

For Model D’s orbit, however, it seems un¬ 
likely that the streamers could also be ob¬ 
served as moving groups in proper motion 
studies. The square-degree surveyed at 40° 
contains particles at ~ 37 kpc with tangential 
velocities of ~ 370 km/s, which corresponds 
to a proper motion of 2.2 X 10 -3 arcsecs/year, 
well below the current observational limits for 
such surveys. 

3.7. Moving Groups following Disrup¬ 
tion 

In Figure 10 we consider the persistence of 
moving groups in the halo, following Model D 
beyond its complete destruction on its sixth 
pericentric passage. The model maintains a 
“core” of unbound stars at densities higher 
than the tidal streamers for more than 1 Gyr 
after this passage. This “core” could not be 
confused with a dwarf spheroidal galaxy, how¬ 


ever, since its surface density quickly falls be¬ 
low 0.1M@/pc 2 , an order of magnitude smaller 
than that of the Galaxy’s satellites (Mateo et 
al. 1993). The persistence of this “core” does 
not contradict estimates of the timescale for 
the dwarf to disperse following its disruption 
(~ r 0 /cr ~ (500 pc)/(10 km/s) ~ 10 8 years). 
At 1 Gyr after disruption the dwarf’s surface 
density has already fallen by a factor of 10 
because of the free streaming of the higher 
velocity stars away from the region. The re¬ 
maining “core” consists of the lowest velocity 
stars whose dispersion is much less than 10 
km/s and, hence, whose dispersion timescale 
is much longer than 10 8 years. 

Panels (a) and (b) of Figure 10 show the 
velocity distribution of stars in the “core” in 
two different magnitude bins, 0.5 Gyrs after 
pericenter. Panel (c) is a similar plot, 1 Gyr 
after pericenter. All three suggest that such 
a group could be observed. Panel (d) is the 
velocity distribution of stars 10° further along 
the orbit than those in panel (c), and here the 
moving group is not so obvious. Note, how¬ 
ever, that these stars are more than 160 kpc 
away from the center of the Galaxy. Presum¬ 
ably it would be less challenging to detect a 
disrupting dwarf on a more tightly bound or¬ 
bit. 

3.8. Carbon Stars and Globular Clus¬ 
ters 

The Fornax dSph contains ~65 carbon 
stars (Richer & Westerlund 1983) and 5 glob¬ 
ular clusters (Hodge 1988). The disruption of 
such a system might leave some signatures of 
the event as substructure in the halo carbon 
star and globular cluster populations. Indeed, 
IGI report the existence of four globular clus¬ 
ters in the vicinity of Sgr itself. 

In the previous section, the square-degree 
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examined in Figure 10(a) contains ~5% of the 
initial mass of the model. For a Fornax-like 
population, such a region would contain ap¬ 
proximately 3 high-velocity carbon stars, giv¬ 
ing a surface density significantly higher than 
the average (0.02 deg -2 ) estimated to a depth 
of M v = 18 for faint, high-latitude carbon 
stars by Green, Margon, Anderson & Cook 
(1994). 

Clearly, the strategy of identifying debris 
by its over-density in a small region of phase- 
space would not work for globular clusters 
from disrupted dwarfs. However, the satellite 
galaxies of the Milky Way tend to lie in the 
great circles defined by the Magellanic Stream 
(the Magellanic Plane) and the Fornax-Leo- 
Sculptor stream, as noted first by Lynden- 
Bell (1976, 1982). Globular clusters have 
subsequently been associated with both these 
planes (Lin & Richer 1992; Majewski 1994). 
These great circles present obvious targets to 
search for moving groups and their discov¬ 
ery would be further proof that the alignment 
of halo objects is not mere coincidence, but 
rather the signature of tidal interactions and 
accretions of satellites by the Milky Way. 

4. Discussion 

Five main conclusions follow from our work: 

1. For pericenters of ~ 13 kpc the con¬ 
tribution of the halo to the disruption 
of a dwarf spheroidal of mass less than 
10 8 M@ is at least that of the disk. 

2. The dwarf galaxy in Sagittarius is nei¬ 
ther necessarily on its first nor on its 
last Galactic passage, despite tidal ra¬ 
dius arguments. 

3. From existing observations it is not pos¬ 
sible to distinguish amongst a wide range 


of orbital parameters and initial dwarf 
models. It is also not possible to de¬ 
termine the direction of motion of Sgr 
from its distorted state and, so, we can¬ 
not say if we are currently observing it 
immediately prior to or following disk 
passage. 

4. Debris from the dwarf’s most recent 
passage should be observable as moving 
groups in the halo. 

5. It may be possible to infer whether or 
not similar events have occurred within 
the past 10 9 years from observations of 
moving groups in the halo. 

Conclusion (1) has implications for mod¬ 
els of Galaxy formation. The existence of the 
thin disk is often used as evidence against sub¬ 
stantial recent accretion of satellites by the 
Galaxy (e.g. Toth & Ostriker 1992; Quinn, 
Hernquist & Fullagar 1993, and references 
therein). In our simulations, if half of the 
initial binding energy of a 10 8 M@ satellite 
is converted entirely into random energy in 
a 10 n M@ disk, this would increase the z- 
velocity dispersion ( a z ) in the disk by only 
~ 0.25 km/s. Moreover, the important con¬ 
tribution of the halo to tidal disruption im¬ 
plies that dwarfs could easily be destroyed on 
orbits with larger pericenters by halo shock¬ 
ing alone, with negligible effect on the disk. 

These findings are consistent with earlier 
studies which stress the fragility of disks (Quinn 
and Goodman 1986; Quinn, Hernquist & Ful¬ 
lagar 1993; Walker, Mihos & Hernquist 1994). 
In the simulations of Walker et al. (1994), 
for example, the satellite orbit decays quickly 
both radially and into the plane of the disk 
through dynamical friction with the disk and 
the halo. A substantial fraction of both the 
orbital and internal energy of the satellite is 
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lost to the disk and a z can increase signifi¬ 
cantly through a single accretion event. In 
our calculations, the ratio of satellite mass to 
disk mass is an order of magnitude smaller 
than previous studies. Since frictional de¬ 
celeration is proportional to the mass of the 
satellite (M sat ), the timescale for orbital de¬ 
cay (oc 1/M sat ) is much longer (see Binney 
& Tremaine 1987). More of the orbital en¬ 
ergy is carried away by stripped particles and 
less is deposited in the disk. Hence, our esti¬ 
mate of an increase of 0.25 km/s in a z sug¬ 
gests that the Galaxy could absorb 10 9 — 
10 lo M@ in small accretion events and still 
maintain a thin disk. This supports Searle 
and Zinn’s view of the formation of the Milky 
Way from primordial fragments. (For further 
discussion, see e.g. Hernquist [1990b], Hern- 
quist & Quinn [1993], and Majewski [1993b]) 
However, fully self-consistent simulations are 
needed to precisely determine the relative im¬ 
portance of the disk and halo to the orbital 
decay and disruption of small satellites. 

Conclusions (4) and (5) emphasize the value 
of observational kinematic studies. A number 
of moving groups have already been reported 
in the literature (see Preston, Beers & Sheet- 
man 1994; Majewski, Munn & Hawley 1994, 
and references therein). The discovery of 
moving groups in the halo associated with Sgr 
would constrain its orbit, its origin, and its 
history. Observations of such groups in gen¬ 
eral (or limits on their occurrence) would tell 
us something about the frequency of accretion 
events in the recent history of our Galaxy. A 
smooth phase-space distribution argues in fa¬ 
vor of a relatively uneventful past, while a 
non-uniform one points to the more common 
occurrence of dwarf galaxy accretion, much as 
we are witnessing today in Sagittarius. 


Since we do not know the history of the 
Galaxy, we cannot rely on surveys covering 
only a few square arc-minutes (e.g. Bahcall, 
Flynn, Gould & Kirhakos 1994) to tell us 
much beyond the local properties of the region 
observed. If the halo was indeed formed from 
many accretion events we would expect non¬ 
uniformity on square-degree scales because of 
this history and it would only be through a 
survey covering many square-degrees on the 
sky that we could obtain global estimates of 
its content and structure 
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Theory Grant NAGW-2422, NSF Grants 90- 
18526, ASC 93-18185 and AST 91-17388, and 
the Presidential Faculty Fellows program. 
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Model 

M 

^0 

N 

Orbit type 

T* tidal 

Impulse 

A 

10 7 

0.6 

10 5 

1 

0.36 

3.67 

B 

10 8 

0.6 

10 5 

1 

0.77 

0.37 

C 

10 8 

0.6 

10 5 

2 

0.82 

0.26 

D 

10 8 

0.6 

5 x 10 6 

3 

0.77 

0.37 


Table 1: Parameters of simulations. Length scales are in kpc and mass is in M@. N is the number 
of particles used in the simulation and the impulse is expressed in units of the initial binding 
energy of the dwarf galaxy. 


Orbit type 

Pericenter 

Apocenter 

Period 

1 

13.4 

81.5 

1.08 

2 

15.5 

81.5 

1.18 

3 

14.4 

185.4 

2.48 


Table 2: Parameters of orbits. Length is in kpc, time is in Gyrs. 
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Figure Captions 

Figure 1 Orbits in the X-Z plane. The disk 
of the Milky Way lies in the X-Y plane. 

Figure 2 Fractional mass remaining bound 
as a function of time in orbital periods (T) 
for each simulation. The “Fig” labels show 
the times at which the state of a model is 
presented in the following Figs 3-7. 

Figures 3-7 Summary of analysis at times 
indicated in Figure 2. 

(a) Surface density contours. The bold con¬ 
tour represents 1.0M@/pc 2 . Consecutive con¬ 
tours are a factor of y/lO apart. 

(b) Line of sight average velocity (open sym¬ 
bols) and velocity dispersion (closed symbols) 
along the projected major axis. 

(c) Normalized Gauss-Hermite moments (see 
text). 

(d) Line of sight velocity distribution of par¬ 
ticles in the bin at —1° (used in (b) and (c); 
f is the fraction of particles per unit a). The 
solid line is a Gaussian of the same veloc¬ 
ity dispersion. The dotted line includes the 
higher Gauss-Hermite moments. The crosses 
indicate the points from simple binning of the 
data. 

(e) As (f), but for the —9° bin. 

(f) Surface density profile (solid line). For 
comparison, a King model of the given con¬ 
centration (c= !og(rtidal/r core )) is shown by 
the dotted line. 

(g) Point along orbit at which the analysis 
was done. 

Figure 8 Galactocentric radial velocity vs. 
radius for random sample of 2500 particles. 
The solid line is the locus of the orbit in this 
plane. 


Figure 9 Analysis of tidal streamers. N is 
the number of stars per (degree) 2 per 1 km/s 
bin at 40° (panels (a) and (b)) and 75° (pan¬ 
els (c) and (d)) from the dwarf. The dotted 
line shows the background towards the North 
Galactic Pole. The solid line is from the simu¬ 
lation. The different panels are for the stated 
magnitude bins. 

Figure 10 As figure 9, but for the dwarf at 
0.5 Gyrs (panels (a) and (b)) and 1.0 Gyrs 
(panels (c) and (d)) after disruption. 
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